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Abstract 

In the early 2000's, Gourley (2000), Wu et al. (2001), Ashwin et al. (2002) initiated 
the study of the positive wavefronts in the delayed Kolmogorov-Petrovskii-Piskunov- 
Fisher equation 

u t (t, x) = Au(t, x) + u(t, x)(l - u(t -h,x)), u>0, x G R m . (*) 

Since then, this model has become one of the most popular objects in the studies 
of traveling waves for the monostable delayed reaction-diffusion equations. In this 
paper, we give a complete solution to the problem of existence and uniqueness 
of monotone waves in equation (*). We show that each monotone traveling wave 
can be found via an iteration procedure. The proposed approach is based on the 
use of special monotone integral operators (which are different from the usual Wu- 
Zou operator) and appropriate upper and lower solutions associated to them. The 
analysis of the asymptotic expansions of the eventual traveling fronts at infinity is 
another key ingredient of our approach. 

Key words: KPP-Fisher delayed reaction-diffusion equation, heteroclinic 
solutions, monotone positive traveling wave, existence, uniqueness. 
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1 Introduction and main results 

It is well known that the traveling waves theory was initiated in 1937 by Kol- 
mogorov, Petrovskii, Piskunov [20] and Fisher [13] who studied the wavefront 
solutions of the diffusive logistic equation 

u t (t, x) = Au(t, x) + u(t, x)(l - u(t, x)), u>0, x G R m . (1) 

We recall that the classical solution u(x,t) = <fi{v ■ x + ct), \\u\\ = 1, is a 
wavefront (or a traveling front) for (1), if the profile function is positive and 
satisfies <f>(— oo) = 0, 0(+oo) = I. 

The existence of the wavefronts in (1) is equivalent to the presence of positive 
heteroclinic connections in an associated second order non-linear differential 
equation. The phase plane analysis is the natural geometric way to study 
these heteroclinics. The method is conclusive enough to demonstrate that 
(a) for every c > 2, the KPP-Fisher equation has exactly one traveling front 
u(x, t) = <f>(v-x+ct); (b) Eq. (1) does not have any traveling front propagating 
at the velocity c < 2; (c) the profile (f> is necessarily strictly increasing function. 

The stability of traveling fronts in (1) represents another important aspect of 
the topic: however, we do not discuss it here. Further reading and relevant 
information can be found in [6,21,28,36]. 

Eq. (1) can be viewed as a natural extension of the ordinary logistic equation 
u'(t) = u(t)(l — u(t)). An important improvement of this growth model was 
proposed by Hutchinson [18] in 1948 who incorporated the maturation delay 
h > in the following way: 

u'(t)=u(t)(l-u(t-h)), u>0. (2) 

This model is now commonly known as the Hutchinson' s equation. Since then, 
the delayed KPP-Fisher equation or the diffusive Hutchinson's equation 

u t (t, x) = Au(t, x) + u(t, x)(l - u(t -h,x)), u>0, x G R m , (3) 

is considered as a natural prototype of delayed reaction-diffusion equations. It 
has attracted the attention of many authors, see [2,4,11,14,15,17,22,33,35,37]. 
In particular, the existence of traveling fronts connecting the trivial and pos- 
itive steady states in (3) (and its non-local generalizations) was studied in 
[2,4,7,11,16,27,33,35]. Observe that the biological meaning of u is the size of 
an adult population, therefore only non-negative solutions of (3) are of inter- 
est. It is worth to mention that there is another delayed version of Eq. (1) 
derived by Kobayashi [19] from a branching process: 

u t (t, x) = Au(t, x) + u(t - h, x){l - u(t, x)), u>0, x G M m . 
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However, since the right-hand side of this equation is monotone increasing 
with respect to the delayed term, the theory of this equation is fairly different 
(and seems to be simpler) from the theory of (3), see [30,35,38]. 



This paper deals with the problem of existence and uniqueness of monotone 
wavefronts for Eq. (3). The phase plane analysis does not work now because of 
the infinite dimension of phase spaces associated to delay equations. Recently, 
the existence problem was considered by using two different approaches. The 
first method, which was proposed in [35], uses the positivity and monotonicity 
properties of the integral operator 



where (H(fi)(s) = 4>(s)((3 + 1 — 0(s — h)) for some appropriate j3 > 1, and 
e' = e(r 2 — r x ) with T\ < < r 2 satisfying ez 2 — z — (5 = 0, and e" 1 / 2 = c > is 
the front velocity. A direct verification shows that the profiles G C(R, R + ) of 
traveling waves are completely determined by the integral equation A<p = <fi. 
Wu and Zou have found a subtle combination of the usual and the Smith and 
Thieme nonstandard orderings on an appropriate profile set T* C C(R, (0, 1)) 
which allowed them (under specific quasimonotonicity conditions) to indicate 
a pair of upper and lower solutions ^ such that _ < A> +1 <f> + < A j (f) + , j = 
0, 1, . . . Then the required traveling front profile is given by = lim A^(j) + . 
More precisely, in [35, Theorem 5.1.5], Wu and Zou established the following 

Proposition 1 For any c > 2, there exists h*(c) > such that if h < h*(c), 
then Eq. (3) has a monotone traveling front with wave speed c. 

The above result was complemented in [33, Remark 5.15] and [27], where it 
was shown that Proposition 1 remains valid if c = 2. It should be observed 
that Wang et al. [33] have also used the method of upper and lower solutions, 
however their lower solution is different from that in [35]. Recently, Ou and 
Wu [26] showed that Proposition 1 can be proved by means of a perturbation 
argument (considering h > as a small parameter). 

The second method was proposed in [11]. It essentially relies on the fact that, in 
a 'good' Banach space, the Frechet derivative of lim e ^ A along a heteroclinic 
solution ip of the limit delay differential equation (2) is a surjective Fredholm 
operator. In consequence, the Lyapunov-Schmidt reduction was used to prove 
the existence of a smooth family of wave solutions in some neighborhood of 
■ip. The following result was proved in [11, Corollary 6.6.]: 

Proposition 2 There exists c* > such that if < h < 1/e then for any 
c > c* , Eq. (3) has a wave solution u(x,t) = 4>{u ■ x + ct), \u\ = 1, satisfying 
0(-oo) = 0, 0(+oo) = 1. 



(A<f>)(t) = - \j e r ^- s \n<j>)(s)ds + J e 



f t +oo 



{t - s) (H<p){s)ds 



(4) 
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We remark that the positivity of this wave was not proved in [11] and the value 
of c* > was not given explicitly. Nevertheless, as it was shown in [12] for the 
case of the Mackey-Glass type equations, the method of [11] may be refined 
to establish the existence of positive wavefronts as well. Moreover, it follows 
from [12] that Proposition 2 is still valid for h G (0,3/2). The recent work [3] 
suggests that the approach of [11] can be also used to prove the uniqueness (up 
to shifts) of the positive traveling solution of (3) for sufficiently fast speeds. 

In this paper, motivated by ideas in [9,35], we give a criterion for the existence 
of positive monotone wavefronts in (3) and prove their uniqueness (modulo 
translation). In order to do this, instead of using operator (4) as it was done 
in all previous works, we work with different integral operators, namely: 

{Ap)(t) = j (e x ^ - e^MsMs - h)ds, (5) 

where e G (0, 0.25) and < A < /i are the roots of ez 2 — z + 1 = 0, and with 

(Bip)(t) = 4 J (s - t)e 2( '- s V(sMs - h)ds (6) 
t 



which can be considered as the limit of A when e — > 0.25. Remarkably, all 
monotone wavefronts (in particular, the wavefronts propagating with the min- 
imal speed c = 2) can be found via a monotone iterative algorithm which uses 
A, B and converges uniformly on R. 

Before stating our main results, let us introduce the critical delay hi = 
0.560771160 . . . This value coincides with the positive root of the equation 



2h 2 exp(l + Vl + Ah 2 - 2h) = 1 + \fl + Ah 2 

and plays a key role in the following result (which is proved in Section 2): 

Lemma 3 Lett G (0,0.25], h > 0. Then the characteristic function tp(z, e) : = 
ez 2 — z — exp(-zh) has exactly two (counting multiplicity) negative zeros Ai < 
A2 < if and only if one of the following conditions holds 

(1) < h < 1/e, 

(2) e > e*{h) and l/e<h< h x . 

Here the continuous e*(h) is defined in parametric form by 

2t 



e*(h(t)) = th(t), hit) = (2t+ v / 4t 2 + 1) exp(-l ), * G [0, 0.445 . . .]. 

1 + v At 2 + 1 

Let us state now the main results of this paper. 
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Theorem 4 Eq. (3) has a positive monotone wavefront u = ^p{y • x + ct), 
\u\ = 1, connecting with 1 if and only if one of the following conditions holds 

(1) < h < 1/e = 0.367879441... and 2 < c < c*(h) : = +oo; 



(2) l/e<h< h = 0.560771160 .. . and2<c<c*(h) :=l/^Je*(h). 

Furthermore, set <j)(s) := ip(cs). Then for some appropriate 0_ (given below 
explicitly), we have that = lim^ +00 A>(j)- (if c> 2), and = lim.,^ +00 B^4>- 
(if c = 2), where the convergence is monotone and uniform on R. Finally, for 
each fixed c ^ c*{h), <f>(t) is the only possible profile (modulo translation) and 
<f>(t), 4>-{t) have the same asymptotic representation 1 — e A2 *(l + o(l)) at +oo. 

Corollary 5 If h > h± = 0.560771160 . . . then the delayed KPP-Fisher equa- 
tion does not have any positive monotone traveling wavefront. 

Next, let us define the continuous function e # (/i) parametrically by 



*/w t + 2 + y/2t + 4 , , , ln(2 + v / 2TT4) , rt , nrt „ 
e # (/i(t)) = ^ , /i(t) = ^ ^ i, t E (-2, -1.806 . . .] (7) 

Set ho '■— 0.5336619208 . . . (see also Lema 8 for its complete definition) and 



c*{h) :-- 



+oo, 



when h E (0,0.5 In 2], 
l/y/e#(h), when h E (0.5 In 2, h ], 
2, when h > ho. 
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Fig. 1. Schematic presentation of the critical speeds and delays. 



Theorem 6 Let u = ip{y ■ x + ct), \v\ — 1, be a positive monotone traveling 
front of Eq. (3). Set <f>(s) := ip(cs). Then, for some appropriate t , positive Kj 
and every small positive a, we have at t = — oo 
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<f)(t + t ) = < 



-K 2 te xt + 0(e ( - 2X -^ t ), whenc = 2, 

e xt - Kxd* + 0(e( 2A - CT )*), when 2 < c < 1.5^2, 



e Xt + 0( e (2A- < r)t) > 

Similarly, at t = +00 



0(t + 1 ) 



wkn c > 1.5v^= 2.121... 

w/ien h < h , c E [2,c*(h)]f] 
when h e (0.5 In 2, /ij] 
and c G (c*(h),c*(h)), 



1 - e A2 * + 0(e( 2A2 +' 7 )*) 
1 - e A2 * + iT 3 e Al *+ 

+0( e (Ai-<r)t) > 

1 - ir 4 te A2t + Oie^ 2 -^), when c = c*(h) and h e (1/e, /ii]. 



(8) 



Theorem 6 suggests the way of approximating the traveling front profile: e.g., 
for c 7^ 2, c*(/i), we can take functions a_(t) := Cie~ A * and a+(t) := 1 — e A2 ' 
and glue them together at some point r. The point r and c\ > have to be 
chosen to assure maximal smoothness of the approximation at r. As we will 
see in Section 3, this idea allows to construct reasonable lower approximations 
to the exact traveling wave. See also Figure 2 below. 

Remark 7 As it was showed by Ablowitz and Zeppetella [1], equation (1) 
has the explicit exact wavefront solution u = (p*(v ■ x + ct), \u\ = 1, with 
c = 5/\/6 = 2.041 . . . and the (scaled) profile 



11 5s \ 2 

- + -tanh(— + s )J , Ms) ■■= <P*(cs). 



If we select so = 0.51n2 ; then 

Ms) = 1 - 2e- 5s/6 ~ 2s » + 0(e~ 5s/4 ) = 1 - e- 5s/6 + 0(e" 5s/4 ), s ->• +00, 

so that (ft* = \im.j^ +OQ A 1 * <fi- in view of Theorem 4 and the uniqueness (up to 
translations) of the traveling front for the non-delayed KPP-Fisher equation. 
Figure 2 (on the left) shows five approximations A-*(f>-,j = 0, 1, 2, 3, 4, and the 
exact solution <f)+, the graphs are ordered as <fi- < A<p- < A 2 (f)- < A 3 4>- < 0*. 
On the right, the four first approximations B^(f)-,j = 0, 1,2,3, of <fi are plotted 
when c = 2, h — 0.56. It should be noted that the limit function <fi and the 
initial approximation </>_ have the same first two terms (1 — exp(A 2 t)) of their 
asymptotic expansions at +00. See Theorem 4 and Sections 3,4- However, as 
the analysis of the Ablowitz- Zeppetella solution shows, these <fi and 0_ may 
have different first terms of their expansions at —00. This partially explains 
a better agreement between the exact solution and their approximations for 
t > t — 0.487. . . on the left picture (the value of r is given in Section 3). 
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Fig. 2. On the left: increasing sequence of approximated waves A 3 '<!>-, j = 0, 1, 2, 3, 4, 
and the Ablowitz-Zeppetella exact solution (e = 0.24 and h = 0). On the right: 
approximations B J (f)-,j = 0, 1, 2, 3 (e = 0.25 and h = 0.56). 

The structure of the remainder of this paper is as follows. In Section 2, the 
characteristic function of the variational equation at the positive steady state 
is analyzed. In the third [the fourth] section, we present a lower [an upper] 
solution. Section 5 contains some comments on the smoothness of upper and 
lower solutions. Theorems 4 and 6 are proved in Sections 6 and 7, respectively. 



2 Characteristic equation at the positive steady state 

In this section, we study the zeros of if)(z, e) := ez 2 — z — exp(-zh), e, h > 0. 
It is straightforward to see that ip always has a unique positive simple zero. 
Since ip"'(z, e) is positive, ip can have at most three (counting multiplicities) 
real zeros, one of them positive and the other two (when they exist) negative. 
Lemma 3 in the introduction provides a criterion for the existence of two 
negative zeros Ai < A 2 < 0. We start by proving this result: 

PROOF. [Lemma 3] Consider the equation — z = exp(-zh). An easy analysis 
shows that (i) this equation has exactly two real simple solutions Z\ < z 2 < 0, 
z 2 > — e, if h G (0, 1/e), (ii) it has one double real root z± — z 2 — — e if 
h— 1/e, and (iii) it does not have any real root if h > 1/e. As a consequence, 

ez 2 — z = exp(-zh) (9) 
has two negative simple solutions if e > and h G (0, 1/e]. 



A similar argument shows that for every h > 1/e there exists e*(h) > 



such that Eq. (9) (a) has two negative simple roots if e > e*(h), (b) has one 
negative double root if e = e*(h), (c) does not have any solution if e < e*(h). 
In particular, e = e*(h),z = Xi(h) = A 2 (/i), solve the system 



ez 2 — z = exp(-zh), 2ez — 1 = —hexp(-zh), 

which yields the parametric representation for e*(h) given in the introduction. 

Finally, a direct graphical analysis of (9) shows that e*(h) is increasing with 
respect to h. Hence, since e*(h) < 0.25, we conclude that h < (e*) _1 (0.25) =: 
h x = 0.560771 ... □ 



Lemma 8 Let X 1 < A 2 < be two negative zeros of ip(z,e) and e G (0,0.25] 
be fixed. Then Ai < 2A 2 if and only if one of the following conditions holds 

(1) < h < 0.5 ln2 = 0.347...; 

(2) e > e*(h) and 0.5 In 2 < h < h := 0.5336619208.... 



PROOF. This lemma can be proved analogously to the previous one, we 
briefly outline the main arguments. First, for each fixed positive e* we may 
find h(e*) > such that A x < 2A 2 if h e (0, h(e*)) and Ai = 2A 2 if h = h(e*). 
In this way, 

e*\\ - A 2 = exp(-A 2 /i(e # )), 4e # A^ - 2A 2 = exp(-2A 2 /i(e # )), 

which yields representation (7). Now, we complete the proof by noting that 
h(e) is continuous and strictly increasing on (0, +oo) and h(0+) = 0.5 In 2, h = 
h{0.25). □ 



Lemma 9 Let Ai < A 2 < be two negative zeros ofip(z,e) and e G (0,0.25] 
be fixed. Then 3?Aj < Ai for every complex root ofip(z,e) = 0. 

PR OOF. S et a := (l+2e-y/i + 4e 2 )/(2e), a := -e~ ah /{Vl + 4e 2 -2e), k : = 
e/(Vl + 4e 2 - 2e). Then a, k > 0, a < 0, and 

rj}{ z + a) = (VI + 4e 2 - 2e)(kz 2 -z-l + ae~ zh ). 

It is easy to see that p(z) := kz 2 — z — 1 + ae~ zh also has two negative and one 
positive root. Since the translation z — > z + a of the complex plain does not 
change the mutual position of zeros of ip, the statement of Lemma 9 follows 
now from [31, Remarks 19,20]. □ 
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3 A lower solution when Ai < A 2 



In this section, we assume either condition (1) or condition (2) of Theorem 
4 holds. In addition, let c G [2,c*(h)) so that A x < A 2 (where Ai := — oo if 
h — 0) and A < fi. Set 

■' -r¥e X(t ~ T \ if t < r, 
1 A I A ~ A2 

r = v ln ^ — T >0 ' M*) 

A 2 A — A 2 



1 



if i > T. 



It is easy to see that 0_ G C^R) nC 2 (M\ {r}) with > 0, tel, and 

e0 / !(t)-0 / _(t) + 0_(t)(l-0_(t-/i)) <0, teR\(r,r + /i]. (10) 
Lemma 10 Inequality (10) holds for all tel. 



PROOF. The case /i = is obvious, so let /i > 0. It suffices to consider 
t G (r, r + /i]. If we take t G (r, r + h], then 



e<^(f) - <t>'_{t) + 0_(t)(l - - h)) = -t\\e Mt + A 2 e A2 '+ 



•V i/i j_ A2 c \(t-r-h)j = _ e x 2 (t-h) e A 2*)(i _|_ 

A — A 2 



(l-e A2 *)(l + 



A, 



A - A 9 



,\{t-r-h) 



) = 



1 _ e Aa(t-h) + A 2 e A(t-T-ft) _ e A 2 t _ c \ 2 t A 2 e A(t-T-ft) 

A — A 2 A — A 2 



1 + 



A 2 



. e A, _ A 6 a 2S _ A _ AA. eA2(s+ , )eAs =; 

A — A 2 A - A 2 A - A 2 (A - A 2 J 2 



where s = t — r — he (—h, 0]. The direct differentiation shows that 
-A 2 A 



As) = 



A-A 9 



A- A, 



>0, 



since e Xs < 1, e A2S > 1, and (1 + 



A2S 



A + A 



2 ^As 



A - A^ 



e As ) > 1, if A + A 2 > 0, 



(l + ^e As )>l + A + A2 



2A 



> 0, if A + A 2 < 0. 



A — A 2 A — A 2 A — A 2 

Finally, we have that p(s) < for all s G [—h, 0] since p'{s) > and 



p(0) = -A 2 (A- A 2 )- 2 e A2h < 0. 



□ 



Remark 11 (A lower solution when Ai = A 2 ) We can not use 0_ as a 

lower solution when c — c*(h), 1/e < h < hi. Indeed, by Theorem 6, in 
this case </>_ converges to the positive steady state faster than the heteroclinic 
solutions. In Section 5, we will present an adequate lower solution for this 
situation. However, it will not be C l -smooth. 
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4 An upper solution when Ai < A 2 

Suppose that Ai < A 2 and set fait) :— 1 — e X2t + e rt for some r e (Ai,A 2 ). 
Recall that Ai := — oo if h — 0. Obviously, ^(r, e) > and 2 (t) £ (0, 1) for 
£ > 0. Next, it is immediate to check that fa : R — >■ R has a unique critical 
point (absolute minimum) to = fair) > 0: 



toir) = 



ln( - r) ~ ln( - A2) , A 2 e^=re-o. 
A 2 - r 



Observe that if h e (0, 1/e), then we can assume that t (r) > h since 

lim toir) = -1/A 2 > l/e> h, 

r-*r\2- 

where the last inequalities were established in the proof of Lemma 3. It is clear 
that the function 

'fait), if t > fair), 
faitoir)), i£t<t (r) 
is C 1 -continuous and increasing on R. Moreover, 0+(i) G C 2 (R\ {toir)}). 



fait) 



Lemma 12 For all r < A 2 sufficiently close to A 2; 0+ satisfies the inequality 
efa'it) - 4>'(t) + 0(0(1 - fat -h))>0, te R. 



PROOF. Step I. First we prove that, for all t > t , the following inequality 
holds: 

imfa)it) := efait) - fa(t) + fait)il - fait - h)) > 0. 

In particular, this implies that (910+) (i) > if t > t + h. For t = t + s, we 
have that 

(9t0 2 )(*) = ^(r, e)e rt - ^(A 2 , e)e A2 ' + (-e A2 * + e ri )(e A2( ^ A) - e r ^) = 
V>(r, e)e r * + (-e A2 * + e ri )(e A2 ('- A ) - e^"^) = 



ipir, e)e rs + ( e AaS + e rs )e rf °(— e^"^ - e r{s ~ h) ) 

A 2 A 2 



3 r(t +s) 



^( r? e ) + (__L e (A 2 -0.5r)s + e 0.5r a ) e rto(iL e -A a ft e (A 2 -0.5r). _ g-rhgO.Sr^ 

A 2 A 2 



V>(r, e) + A 1 (s)e rt M 2 (s) 



It is easy to see that A,-(+oo) = and that Aj has a unique critical point Sj, 
with 

lim Si(r) = — 1/A 2 , lim s 2 (r) = /i — 1/A 2 . 

Therefore, for some small 5 > and for all r close to A 2 , the function 
Ai(s)e rto A 2 (s) is strictly increasing to on the interval [h — 1/A 2 + 5, +oo) 
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and it is strictly decreasing on [0, — I/A2 — 6]. This means that if (9102) (t) > 
for all t e [t -l/\ 2 -6,t + h- 1/A 2 + 5] then (9?0 2 )(t) > for t > t . In 
order to prove the former, consider the expression 



e ~rt 



(e$(t) - 0' 2 (t) + 2 (t)(l - Mt - &))) = 



r - A 2 

Jl fi A 2s _|_ R rs') R rtoC T »L R A2(s-h) _ e r(s-h)} 

■= r e (r, s). 



^(r,e)e rs + (~^e A2S + e rs )e r *°(^e A2 ( s - /l ) - e r(s_fc) ) 



r - A 2 

Since r e (r, s) is analytical on some open neighborhood Q C M 2 of the compact 
segment {A 2 } x [— 1/A 2 — 5, h — 1/A 2 + 5] C M 2 , we find that, for every fixed 
e > 0, 

lim r e (r,s) = ^'(A 2 ,e)e A2S < 
uniformly on [— 1/A 2 — 5, h — I/A2 + 5]. As a consequence, we obtain that 

e0 2 '(t)-0 2 (t) + 2 (t)(l-02(t-/i)) > 0, t e [t -i/\ 2 -8,to + h-l/X 2 + S\. 

Step II. Now, we are ready to prove that (9T(/> + )(t) > 0, t G [t ,to + h]. Indeed, 
since ^2(^0) < 02 (i — ^) for t G [to^o + h], we have that 

(9^+)(0 = e^(t) - 4>' 2 {t) + 02(0(1 " h(t )) > 

(%(t) - 4>' 2 {t) + 2 (t)(l - fa(t -h))>0, t G [t , to + h}. 

Finally, since the inequality (910+) (t) > 0, t < t , is obvious, the proof of the 
lemma is completed. □ 

Remark 13 (An upper solution when Ai = A2) We can not use + as an 

upper solution when c = c*(h), 1/e < h < hi. Moreover, in this case it is 
not difficult to show that + satisfies inequality (10) for all r < A 2 sufficiently 
close to A2 and for large positive t. 



5 Some comments on upper and lower solutions 

5.1 Non-smooth solutions 



The problem of smoothness of the lower (upper) solutions is an interesting 
and important aspect of the topic, see [5,23]. As we have seen in the previ- 
ous sections, C 1 — smoothness condition can be rather restrictive even when a 
simple nonlinearity (the birth function) is considered. The above mentioned 
works [23] show that continuous and piece- wise C 1 — continuous lower (upper) 
solutions 0± still can be used if some sign conditions are fulfilled at the points 
of discontinuity of 0'_j_. Moreover, as we prove it below even discontinuous 
functions <p± can be also used. We start with a simple result of the theory of 
impulsive systems [29] which can be viewed as a version of the Perron theorem 
for piece- wise continuous solutions, cf. [5]. 



11 



Lemma 14 Let (|):R->Kfc« bounded classical solution of the second order 
impulsive equation 



r + atf + b*p = f(t), A^lt, = aj, A^'| t . = ft, 

where {tj} is a finite increasing sequence, f : R — > R is bounded and continu- 
ous at every t ^ tj and the operator A is defined by Aw\ t . := w(tj+) —w(tj—). 
Assume that z 2 + az + b = has two positive roots < A < /i. Then 



+oo 

if A < fj, we have that ^(t) = J (e A( *- s) - e^- s) ) f{s)ds (11) 

+— L- 2 [(Ae^*-*^ - fie^-^) a 3 + (e A ^) - e**"*')) ft] , t ^ t j; 



if A = fi = -0.5a we have that ip(t) = J (s - t) e -°- 5a ^- s) f(s)ds + 

t 



t<tj 



PROOF. See [29, Theorem 87]. □ 

Next, the corollary below shows that our lower solution is an upper solution 
in the sense of Wu and Zou [35] : 

Corollary 15 Assume that ip : R — > R is bounded and such that the deriva- 
tives ip',ip" : R \ {tj} — > R exist and are bounded. Suppose also that ip is a 
classical solution of the impulsive inequality 

r + atf + bij< fit), A^| t . = a 3 , M>\. = ft. 

If oij > 0, ft < 0, then 

+oo 



i>{t) < — — ^ J (e A( '~ s) - e^- s) ) f(s)ds, when A < fi, 
^ t 

ip(t) < J (s - t)e-°- 5a ^- s) f(s)ds, when A = fj, = - 



0.5a. 



PROOF. Suppose that A < /i, the case A = /x is similar. Clearly, q(t) : = 
/(f) - (ip"(t) + ai)'{t) + bip(t)) > and Ae"^) < fie x ^\ e^-^ > e^*"^ 
for t < tj. Thus the desired inequality follows from (11). □ 
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5.2 A lower solution when Ai = A2, h e (1/e, /ii] 



In Section 3, a lower C 1 — solution was presented for the case when Ai < A2. 
However, to apply our iterative procedure in the critical case Ai = A2, we also 
need to construct a lower solution for the corresponding range of parameters. 
It is worth to mention that our approach does not require any upper solution 
once a lower solution is found and the existence of the heteroclinic is proved, 
see Corollary 26. Here, we provide a continuous and piece-wise analytic lower 
solution <f>-(t) if Ai = A2. Our solution has a unique singular point r' where 
A0_| r / = 0, A0'_| r / > 0. This shows that, in general, the sign conditions of 
Corollary 15 need not to be satisfied. 

Take some positive A > (e~ X2h — l)/h and let t' be the positive root of the 
equation At + 1 = e~ x ' 2t . It is easy to see that r' > h. Consider the piece- wise 
smooth function 0_ : R — >■ [0, 1) defined by 



-(*) = 



0, if t < t', 

1 - (At + l)e A2 *, if t > t'. 



12) 



Proposition 16 The inequality (/C0_)(t) > <j>-(t) holds for all tel. 

PROOF. Below, we are assuming that h 7^ hi so that A < // and /C = A; 
however, a similar argument works also in the case h = hi (when /C = B). It 
suffices to prove that (A(f>-)(t) > 4>-(t) for t > t'. Let C 2 — smooth function 
ip be defined by 



1 - (At + l)e A2 *, if t > t' - h, 
B(t), if0<t<r'-h, 
0, if t < 0, 



for some appropriate continuous decreasing B(t). Set 

C(t) := ef{t) - i//(t) + fl;(t)0- - ^ - h)). 

It is easy to check that ( G C(R, R) is bounded on R and ((t) < for all 
t > t'. But then, for all t > r', we have that 



e(/i - A) 



(< 



X(t-s) _ p n(t-s) 



c(s)d S < (Am) < 



+00 

— — J (e x{t ^ - e^*- s ))0_(s)0_(s - h)ds = (A(f)-)(t). □ 



e(/i - A) 
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5. 3 Ordering the upper and lower solutions 

Finally, we show that the condition of the correct ordering 0_ < <p + is not 
at all restrictive provided that solutions 0± are monotone and satisfy some 
natural asymptotic relations. 

Lemma 17 Assume that functions (j>± : R — >■ [0, 1), j = 1,2, are increasing 
and, for some fixed k G {0, 1} ; the following holds 

lim Mt'K Xt = <x±, lim (1 - ± (t))r fc e- A2i = (3±, k , 

t— >— oo t— >+oo 

where (3± t k > 0, and a_ G [0, +oo) ; a + G (0, +oo]. Then there exists a real 
number a such that <f>-(t) < <f>+(t + a) for all tel. 

PROOF. It is clear that 0_(— oo) = and 0±(+oo) = 1. Let <To be sufficiently 
large to satisfy (5 + ^e X2(T0 < P-,k, a-e~ Xa " < a + . Then there exist ti,t 2 such 
that t\ < t 2 and — a ) < (f)+(t), t G X :— (— oo,£i] U [t 2 ,+oo). Now, set 
a — <7o + (t 2 — ti). Since both functions are increasing, we have 

<p-(t-a) < <f>-(t-cr ) < (j>+(t), tel, 
0_(t-<7)<0+(t-(t 2 -ti))<0+(t), te[hM]- □ 



6 Proof of Theorem 4 

6.1. Necessity. Let u(t,x) = ((ct + v • x) be a positive bounded monotone 
solution of the delayed KPP-Fisher equation. Then ip(t) = ((ct) satisfies 

e<p"(t)-<p'(t) + ip(t)(l-ip(t-h)) = 0, teR, (13) 

t 

e<p'(t) = €(p'(0) - v?(0) + (p(t) + J <p(s)(l - cp(s - h))ds. 



The latter relation implies that (p(±oo) G {0, 1} since otherwise <//(±oo) = oo. 
Hence ip : R ->■ (0, 1). Let G C 2 (1R, (0, 1)) be an arbitrary solution of (13). 
Suppose for a moment that 4>'(to) = 0. Then necessarily <f)"(t ) < so that to is 
the unique critical point (absolute maximum) of <fi. But then <j>'(s) < for s > 
to, so that (f)"(s) < 0, s > t , which yields the contradiction 0(+oo) = — oo. 
In consequence, either (f)'(s) > or (f)'(s) < for all s G R. But as we have 
seen, <f>'(s) < implies 0(+oo) = — oo, a contradiction. Hence, any solution 
G C 2 (R, (0, 1)) of (13) satisfies <p'(t) > 0, 0(-oo) = 0, 0(+oo) = 1. 

Lemma 18 If (p G C 2 (R, (0, 1)) satisfies (13), then e G (0,0.25]. 
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( 1 \ 


f o o\ 






\-l/e l/ej 


^(t-hyeOj 



PROOF. Suppose for a moment that e > 0.25. Then the characteristic equa- 
tion eA 2 — A + 1 = associated with the trivial steady state of (13) has two 
simple complex conjugate roots u± = (2e) _1 (l ± iy/4e — 1). 

Since G C 2 (K, (0, 1)) is a solution of (13), it holds that <f>'(t) > 0, t G K, 
(f>(— oo) = 0. Set = ((f>(t), <f)'(t)) T , it is easy to check that z(t) satisfies the 
following asymptotically autonomous linear differential equation 



z '(t) = (A + R(t))z(t), t G 



Since R(—oo) = 0, /° M \R'(t)\dt = 4>{—h) and the eigenvalues u± of A are 
complex conjugate, we can apply the Levinson theorem [10, Theorem 1.8.3] 
to obtain the following asymptotic formulas at t — — oo: 

<t>{t) = (a + o(l))e*/( 2e ) cos(t v / 4^T(l + o(l)) + b + o(l)), 
0'(t) = ( c + (l))e^ sin(t v / 4^T(l + o(l)) + d + o(l)), 

where a 2 + c 2 7^ 0. But this means that either <f>(t) or (f>'(t) is oscillating around 
zero, a contradiction. □ 

Lemma 19 If h > hi or h G (l/e,/ii] and c > c*(/i) i/ien Eq. (13) does not 
have any solution G C 2 (1R, (0, 1)). 

PROOF. On the contrary, let us assume that Eq. (13) has a solution G 
C 2 (1R, (0, 1)). Then Lemma 18 implies that e G (0,0.25] and therefore the 
assumptions of this lemma imply that ip(z,e) does not have negative zeros. 
Following the approach in [32], we will show that this will force <f>(t) to oscillate 
about the positive equilibrium. For the convenience of the reader, the proof is 
divided in several steps. 

Claim I: y(t) := 1 — <f>(t) > has at least exponential decay as t — >■ +00. First, 
observe that 

ey"(t)-y'(t) = <l>{t)y(t-h), teR. (14) 

Therefore, with 7 := (f)(t ), which is close to 1, and g(t) := (f)(t)y(t — h) — 
(f>{to)y(t), we obtain that 

ey"(t)-y'(t)- 1 y(t)-g(t)=Q J teR. 

Note that g(t) > for all sufficiently large t. Since y(t),g(t) are bounded on 
R, it holds that 

/ t +OO N, 

^ = -^37) / e l{t - s) g(s)ds+ J e m ^g(s)ds 



t 
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where / < and < m are roots of ez 2 — z — 7 = 0. The latter representation 
of y(t) implies that there exists T such that 



+00 



y'(t) - ly(t) = ~ J e m ^g(s)ds < 0, t > T . (15) 

t 

Hence, (y(t) exp(-lt))' < 0, t > T , and therefore 

y(t) < y(s)e l ^ s \ t>s>T , g(t) = 0(e lt ), t +00. (16) 

It is easy to see that these estimates are valid for every negative / > (2e) _1 (l — 
VI + 4e). Finally, (15), (16) imply that y'(t) = 0(e lt ), t +00. 

Claim II: y(t) := 1 — (p(t) > is not superexponentially small as t — )■ +00. 
We already have proved that y(t) is strictly decreasing and positive on R. 
Since the right hand side of Eq. (14) is positive and integrable on R + , and 
since y(t) is a bounded solution of (14) satisfying y(+oo) = 0, we find that 

+00 

y(t) = / (1 - e^/^sMs - h)ds. (17) 

t 

As a consequence, there exists 7\ such that 

t t 
y(t) > 0.5(1 - e~°- 5h/e ) J y(s)ds := £ J y(s)ds, t>T x -h. 



t-0.5h t-0.5h 



Now, since y(t) > for all t, we can find positive C, p such that y(s) > Ce~ ps 
for all s E [7\ — /i, T x ] . We can assume that p is large enough to satisfy the 
inequality £(e°- 5ph -l) > p. Then we claim that y(s) > Ce~ ps for all s>Ti-h. 
Conversely, suppose that t' > Ti is the leftmost point where y(t') = Ce~ pt . 
Then we get a contradiction: 

t' t' 

r- r- p 0.5ph _ 1 

y(f) > Z / y(s)ds >Ci j e~ ps ds = Ce~ pt 'i > Ce~ pt ' . 

t'-0.5h t'-0.5h ^ 

Claim III: y(t) > can not hold when ip(z, e) does not have any zero in 
(—oo,0). Observe that y(t) = 1 — <j){t) satisfies 

ey"(t) - y'(t) - (1 - y(t))y(t - h) = 0, t G R, 

where in virtue of Claim I, it holds that (y(t),y'(t)) = 0{lt) at t — +00. 
Then [25, Proposition 7.2] implies that there exists 7 < / such that y(t) = 
v(t) + 0(exp(7t)), t — > +00, where v is a non empty (due to Claim II) finite 
sum of eigensolutions of the limiting equation 

ey"{t) - y\t) - y(t - h) = 0, t e R, 
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associated to the eigenvalues Xj G F — {7 < 9?Aj < /}. Now, since the set 
F does not contain any real eigenvalue by our assumption, we conclude that 
y(t) should be oscillating on R + , a contradiction. □ 

6.2. Sufficiency. Suppose that e G (0,0.25] and let < A < /i be the roots 
of the equation ez 2 — z + 1 = 0. In Lemmas 20-23 below, /C stands either for 
A or B (defined by (5), (6)). 

Lemma 20 If<p, ip G C{R, (0, 1)) and <f>(t) < i/>(t) for all t G R, then /C0, Ktj} G 
C(R, (0, 1)) and (K<j))(t) < (K,ip)(t), tel. Moreover, if increasing then 
/C0 is also increasing. 

PROOF. The proof is straightforward. □ 

Lemma 21 Lei e G (0,0.25]. 7/0+ G C 1 (R, (0, 1)) satisfies the inequality 

e<f)"(t) - </>'(t) + 0(t)(l - <f>(t -h))>0 

for all t G R' := R \ {T 1: . . . ,T m } and + (t),0' + (t) are bounded on R' , then 
(/C0+)(t) < 0+(O /or a// i G R. 

PROOF. If w(7i) := and 

:= ect>'[{t) - cf>' + (t) + + (O(i - + (t - /i)), teR' = R\{r 1 ,...,r m } 

then u(t) > for all t G R', u;(£) is bounded on R' and 

60+ (t) - 0' + (*) + + (*) = Wl (t) , £ G R , 

where Wi(£) := w(£) + + (£)0 + (£ — /i) is bounded on R'. Let now e G (0, 0.25). 
By Lemma 14, we obtain that 

+00 

0+(£) = / (e A (*- s ) -e" (t - s) Ws)ds = 

e(/i - A) 7 

+00 

+ ^73^) / (e A(t " a) " e^ s) ) u(s)ds > (Acf> + )(t). 

The case e = 0.25 (which corresponds to K — B) is completely analogous to 
the previous one. □ 

The proof of the next lemma is similar to that of Lemma 21: 

Lemma 22 Let e G (0,0.25]. 7/0_ G C^R, (0, 1)) satisfies the inequality 

e0"(£) - 0'(£) + 0(t)(l - 0(£ - /i)) < 

for allt G R\ {Ti, . . . , T m } and 0" (£), (f>'_(t) are bounded on R \ {Ti, . . . , T m } 7 
£/ien (/C0_)(£) > 0_(£) /or a// £ G R. 
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Set (f)f +1 := (K,(j)f), j > 0, 0o := ± , and let the increasing functions 0^ < 0+ 
be as in Lemmas 21, 22. Then 

0_ < 07 <...<$-<$+<... 07 ...<... <f>t < 0+, 

where $±(t) = lim^^ 0^(i) pointwise and (fif are increasing (by Lemma 20). 

Lemma 23 $± are wavefronts and $±(t) = lim^oo 0^(t) uniformly on R. 

PROOF. Applying the Lebesgue's dominated convergence theorem to 0~ +1 := 
/C0~, we obtain that <&_(£) = (/C$_)(t). Differentiating this equation twice 
with respect to t , we deduce that $_ : R — >■ (0, 1) is a C 2 -solution of (13) 
(and thus $'_(£) > 0). As a consequence of the Dini's theorem, we have that 
$_(t) = linij^oo (j)j(t) uniformly on compact sets. Since 0~ are asymptot- 
ically constant and increasing, this convergence is uniform on R. The proof 
for $ + is similar. □ 

Corollary 24 Eq. (3) has a monotone wavefront u(x,t) = ((x-v + ct), \u\ = 
1, connecting with 1 if one of the following conditions holds 

(1) 0<h< 1/e and 2 < c; 

^) l/e<h< h and2<c< c*(h). 

PROOF. It is an immediate consequence of Lemmas 10, 12, 17, 21-23. □ 

If c = c*(h), the reasoning of the last proof does not apply because of the lack 
of explicit upper solutions. Below, we follow an idea from [32, Section 6]: 

Lemma 25 Eq. (3) has a positive monotone wavefront u(x,t) = ((x- is + ct), 
\u\ = 1, connecting with 1 if 1/e < h < hi and c = c*(h). 

PROOF. Case I. Fix some h € (1/e, hi) and e = e*(h). Then there exists a 
decreasing sequence ej ! e*(/i) such that Eq. (13) has at least one monotone 
positive heteroclinic solution (f>j(t) normalized by 0j(O) = 0.5. It is clear that 
4>j(t) = (A<j>j)(t). Moreover, each yj(t) :— 1 — <j>j(t) > solves (17) so that 

+oo 

|0^.(t)| = I J e (*- s )A0( s )(i - 0^.(5 - h))ds\ < 1, tel. 
t 

Thus, by the Ascoli-Arzela theorem combined with the diagonal method, {4>j} 
has a subsequence {<fij k } converging (uniformly on compact subsets of R) to 
some continuous non-decreasing non-negative function 0*, 0*(O) = 0.5. Apply- 
ing the Lebesgue's dominated convergence theorem to 4>j k {t) = (A<f)j k )(t), we 
find that 0* is also a fixed point of A. Hence, 0* : R — > [0, 1] is a monotone so- 
lution of Eq. (13) considered with e = e*(h). Since 0*(O) = 0.5, 0* : R ->■ (0, 1) 
is actually a monotone wavefront. 
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Case II. Finally, let e = 0.25 and h — hi. This case can be handled exactly in 
the same way as Case I if we keep e = 0.25 fixed, replace A with £>, and take 
some increasing sequence hj t hi instead of ej I e*(/i). □ 

Corollary 26 Assume that c = c*(h), 1/e < h < hi, and let 0_ be as in 
(12). If A is sufficiently large, then 

0- < < • • • < <Pj ■ ■ ■ < $ = 
where $ is a wavefront and $(£) = lim^oo <fij(t) uniformly on R. 

PROOF. If c = c*(/i) we will take the heteroclinic solution whose exis- 
tence was established in Lemma 25 as an upper solution. Due to (8), we can 
assume that 

(3+ i := lim (1 - ^_(t))*- 1 e -A2 * > 0. 

t— >-\-OG 

Next, let <j>_ be defined by (12). Since a_ := lim^.^ (f)_{t)e~ xt = and 

:= lim (1 - - 4))i _1 e- A2 * = Ae"^ > fi +tU 

for sufficiently large A, Lemma 17 implies that (f>-(t) < + cr), i G R, for 
some a. Finally, it suffices to take := $_(£ + cr) and repeat the proof of 

Lemma 23. □ 

6.3. Uniqueness. Our method of proof follows a nice idea due to Diekmann 
and Kaper, see [9, Theorem 6.4]. Suppose that c ^ c*(h) and let 0i,02 be 
two different (modulo translation) profiles of wavefronts propagating at the 
same speed c. Due to Theorem 6, we may assume that 0i,02 have the same 
asymptotic representation <pj(t) = 1 — e A2 '(l + o(l)) at +oo. Moreover, <pj = 
/C0 j ,where/C = „4ifc>2and/C = i3ifc = 2. Set u(t) := \<f> 2 (t) -0i(t)|e- A2 *. 
Then u(±oo) = 0, u(t) > 0, t € R, and w(r) = max seK w(s) := \cj\ > for 
some r. From the identity 2 — <f>i = /C0 2 — /C0i, we deduce that 



w(r) < r- / (e A ( T ~ s ) - e^ T ^)(w(s)e A2S + w(s - h)e x ^ s -^)ds < 



7° e / (e A ( T - s ) - e^ r - s ))(e A2S + e^^ds = |w| = w(r), if c> 2; 
e(/x — X) J 

+oo 

w(r) < 4e- A2T y (s - r)e A(r - s) (w(s)e A2S + w(s - /i)e A2(s - fe) )ds < 

T 

+ 00 

4|w| e- A2T J (s- r)e x(T - s \e X2S + e X2{s ~ h) )ds = \u\ = w(r), if c = 2, 

r 

which is impossible. Hence, \uu\o = and the proof is complete. □ 
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7 Proof of Theorem 6 



First, using the bilateral Laplace transform (Cy)(z) := j R e sz y(s)ds (see e.g. 
[34]), we extend [25, Proposition 7.1] (see also [3, Lemma 4.1] and [32, Lemma 
22]) for the case J = R. 

Lemma 27 Set x(z) := z 2 + az + f} + pe~ zh and let y G C 2 (R, R) satisfy 
y"(t) + oV(0 + Py(t) + pj/(t - h) = f(t), t G R, (18) 



where a, j3,p, h G R and 

C^e"- 8 *), as t ->■ +oo, 



y(0 = 



0(e bt ), as t ->■ -oo; 



/(*) = 



0(e~ c *), as t ->■ +oo, 
0(e rf ), as t — >■ — oo, 



(19) 



/or some non-negative b < c, B < C,b + B > 0. Tnen, /or eacn sufficiently 
small a > 0, it holds that 



y(t) 



w+(t) + e- {c - a)t o(l), as t ->■ +oo, 
+ e (c " a) 'o(l), as i ->• -oo, 



where 



w±(t) = ± £ Res 



2 = A,' 



A,eF± 



/ e-"f( 8 )ds 



zs a finite sum of eigensolutions of equation (18) associated to the eigenvalues 
\ j eF + = {-C + a <$t\i< -B} and Xj G F_ = {& < KA; < c - o}. 



PROOF. We will divide our proof into several parts. 

Step I. We claim that there exist non- negative B', b' such that B' < B,b' < b, 
B' + 6' > and 

0(te- B>t ), as t ->■ +oo, 
y'(t),y"(t) = { (20) 

0(te bt ), ast^-oo. 



We will distinguish two cases: 

Case A . Suppose that a — 0. Then clearly y"{t) = 0(e~ Bt ) at t — — oo, 
is bounded on R and therefore y'(t) is uniformly continuous on R. Since 
B + b > then either w(+oo) = 0, limsup^^^ \y(s)\ < oo or y(—oo) = 
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0, limsup s ^ +OC) \y(s)\ < oo. Suppose, for example that B > (hence y(+oo) = 
0), the other case being similar. Then, applying the Barbalat lemma, see e.g. 
[35], we find that y'(+oo) = 0. This implies that y'(t) = — y"(s)ds = 
0{e- Bt ) at t = +oo. Thus we may set B' = B. Now, y'(t) = j/'(°)+/o v"( s ) ds = 
0(t) at t = — oo so that we can choose b' = 0. 

Case B . Let now a^O. For example, suppose that a > (the case a < is 
similar). Then, for some £, 

t 

y'(t) = ie- at + J e-^- s \f{s) - (3y(s) - py(s - h)}ds. 

— oo 

In fact, since the second term of the above formula is bounded on R and we 
can not have y'(-oo) = ±oo (due to the boundedness of y(t)), we obtain that 
f = 0. But then y'(t) = 0(e bt ), t ->■ -oo and y'(t) = 0(te~ min ^ B ^), t ->■ 
+oo. Note that b' + B' — min{a + 6, 5 + 6} > 0. Finally, (18) assures that 
(20) is also valid for y"(t). 

Step II. Applying the bilateral Laplace transform C to (18), we obtain that 
x(z)y(z) = f(z), where y = Cy, f = Cf and —B' < $lz < b'. Moreover, from 
the growth restrictions (19), we conclude that y is analytic in — B < $tz < b 
while / is analytic in — C < ^Stz < c. As a consequence, H(z) = f(z)/x( z ) is 
analytic in —B < $lz < b and meromorphic in — C < ^Stz < c. Observe that 
H(z) = 0(z~ 2 ), z — > oo, for each fixed strip U(si,s 2 ) = {si < < s 2 }, 
—C < s 1 < s 2 < c. Now, let a > be such that the vertical strips c — 2a < 
$tz < c and — C < $tz < —C + 2a do not contain any zero of x( z )- By the 
inversion formula [34, Theorem 5a], for each 5 G (—B,b), we obtain that 

y(t) = 2~ / e zt y(z)dz = — f e zt H(z)dz = w±(t) + u±(t), t G R, 

<5— ioo <5— ioo 

, ~ v =p(c— cr)+ioo 

where w±(t) = ± £ Res 2=Aj ^^, u±(t) = — f e zt H(z)dz. 

X[Z) Z7r \(c-a)-ioo 

The above sum is finite, since x( z ) nas a finite set of the zeros in F±. Now, 
for a(s) = H(=f(c — a) + is), we obtain that 

u±(f) = 2 ^ ^ y e^a(s)rfs L t G R. 

Next, since a G £i(R), we have, by the Riemann-Lebesgue lemma, that 

lim / e^aAs^s = 0. 
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Thus we get u±(t) = e T<c ~°")*o(l) at t — oo, and the proof is completed. □ 
Now we can prove Theorem 6: 

PROOF. [Theorem 6] Case I: asymptotics at t — +00. It follows from (16) 
that y(t) = 1 — 4>(t) satisfies y(t) = 0(e lt ), t — > +00, for every negative 
I > (2e)- 1 (l - v / T+4e)- Moreover, /(*) := -y(t)y(t - h) = 0(e 2H ),t^ +00, 
y(t) = 0(1), t ->■ -00 and 

ey"(t) - y'(t) - y(t - h) = -y(t)y(t - h), t e R. 

Therefore Lemma 27 implies that, for every small a > 0, 

y(t) = £ Res, =A ,^&i + e^o(l), t +00. 

2l+a<St\i<l X\ Z ) 



Now, observe that (2e) _1 (l - y/1 + 4e) > A 2 so that either A 2 G (2/ + cr, Z) 
or A 2 < 21. In the latter case, we obtain y(t) = e^ 2l+CT ^ t o(l) , i — >■ +00, which 
allows to repeat the above procedure till the inclusion A 2 G (2 J 7 + cr, 2- 7 " 1 /) is 
reached for some integer j. In this way, assuming that Ai < A 2 , for each small 
a > 0, we find that 

y (t) = rje X2t + Oie^t), where rj := /r^'M^M* ~ *0<k > . (21) 

-X (A 2 ) 



Now, if c = c*(/i) (i.e. A x = A 2 ), we obtain analogously that 
y(t + t ) = ite X2t + 0(e( A2 "^), t -> +00, 
for some appropriate t an d £ > 0. 

Suppose now that ft, G (0,/ioL c < c#(h). Then Lemmas 8, 9 imply that 
3?Aj < Ai < 2A 2 . This means that formula (21) can be improved as follows: 

y (t) = r]e X2t + 0(e {2X2+ ^), t ->• +00. 

Finally, if /i G (0.51n2,/i ] and c G (c # (/i), c*(/i)), it holds that 2A 2 < A x < A 2 . 
Then 



y(t) = £ Res 2=Aj ^M + e (^+^o(l) = 

2A 2 +a<JiA : ,<A2 K\ Z > 

^ + 9e xn + e (M-«)t o{ll where 9 ;= /Re- Aia y(a)y(a-fe)da < Q 

-X (AiJ 
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Case II: asymptotics at t — — oo. This case is much easier to analyze since the 
characteristic polynomial ez 2 — z + 1 of the variational equation 



along the trivial equilibrium of (13) has only two real zeros < A < //. It is 
easy to check that 2A < /j, if and only if c > 1.5a/2 = 2.121 

Since <f>(— oo) = and equation (22) is exponentially unstable on R_, we 
conclude that the perturbed equation 



is also exponentially unstable on R_ (e.g. see [8]). As a consequence, <f)(t) = 
0(e mt ), t — > — oo, for some m > 0. Now we can proceed as in Case I, since 



e0"(t)-0'(t) + 0(t) = / 1 (t), 
with fi(t) := (f){t)<j){t -h) = 0(e 2mt ). The details are left to the reader. □ 
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